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^ ■ ABSTRACT 

2 ■ We have observed the polarized hght at 2 in the center of NGC 1068 with iT^T NICMOS 

\ Camera 2. The nucleus is dominated by a bright, unresolved source, polarized at a level of 

■ 6.0± 1.2% with a position angle of 122° ± 1.5°. There are two polarized lobes extending up to 8" 

northeast and southwest of the nucleus. The polarized flux in both lobes is quite clumpy, with 
the maximum polarization occurring in the southwest lobe at a level of 17% when smoothed to 



}_j ■ 0.23" resolution. The perpendiculars to the polarization vectors in these two lobes point back to 



the intense unresolved nuclear source to within one 0.076" Camera 2 pixel, thereby confirming 
that this is the illuminating source of the scattered light and therefore the probable AGN central 
engine. Whereas the polarization of the nucleus is probably caused by dichroic absorption, the 
polarization in the lobes is almost certainly caused by scattering, with very little contribution 
from dichroic absorption. Features in the polarized lobes include a gap at a distance of about 1" 
from the nucleus toward the southwest lobe and a "knot" of emission about 5" northeast of the 
nucleus. Both features had been discussed by ground-based observers, but they are much better 
defined with the high spatial resolution of NICMOS. The northeast knot may be the side of a 
molecular cloud that is facing the nucleus, which cloud may be preventing the expansion of the 
northeast radio lobe at the head of the radio synchrotron-radiation-emitting jet. 

We also report the presence of two ghosts in the Camera 2 polarizers. These had not been 
detected previously (Hines et al. 2000) because they are relatively faint and require observations 
of a source with a large dynamic range. 
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Subject headings: galaxies: individual (NGC 1068) — galaxies: Seyfert — polarization 



1. Introduction 

NGC 1068 is the closest (14.4 Mpc for Ho = 75 km s~^ Mpc~\ 72 pc arcsec"^) Seyfert 2 galaxy; thus 

it has been studied extensively at all wavelengths (sec the summary of the literature by Bland-Hawthorn 
et al. 1997a). The barred spiral galaxy is ~ 7' in diameter with the northern rotation axis apparently 
tipped towards the earth at an angle of about 40° (Schinnerer et al. 2000; Bland-Hawthorn et al. 1997a 
and references therein). A cone of polarized (scattered) optical and ultraviolet (UV) light (e.g., Scarrott et 
al. 1991; Miller, Goodrich, & Mathews 1991; Capetti et al. 1995a, 1995b) and co-spatial [O III] (Pogge 
1988; Bruhweiler et al. 2001) and [N II] (Cecil, Bland, & Tully 1990) emission fan out from the nucleus 
to the northeast at a position angle of ~ 31°). Radio jets emanate from a bright radio "core" along a 
similar position angle (Wilson & Ulvestad 1983), but only a faint [O III] counterpart is observed toward the 
southwest (linger et al. 1992) and very little UV/optical polarized light (Scarrott et al. 1991). The faintness 
of the southwest cone suggests that it lies behind the plane of the galaxy and is extinguished by the dust in 
the plane (e.g., Gallimore et al. 1994; Bland-Hawthorn et al. 1997b), whereas it is thought that the brighter 
northeast cone emerges from the nucleus in front of the galaxy plane (e.g., Kishimoto 1999; Bruhweiler et 
al. 2001). Broad permitted lines including Fe II are observed in the spectrum measured in polarized light, 
confirming that this proto-type Seyfert 2 (narrow-lined) galaxy harbors a Seyfert 1 (broad-lined) nucleus 
(Antonucci & Miller 1985; Antonucci, Hurt, & Miller 1994; Miller, Goodrich, & Mathews 1991; IngUs et 
al. 1995; Alexander et al. 1999) that is hidden from our direct view by a geometrically and optically thick 
dusty region. Light from the nucleus ionizes the clouds of the narrow line region; moreover, ionizing and UV 
photons from the nucleus are scattered by electrons into the more extended, diffuse ionized gas of the galaxy 
disk and halo (Bland-Hawthorn, Sokolowski, & Cecil 1991; Sokolowski, Bland-Hawthorn, & Cecil 1991; Neff 
et al. 1994). The light is also scattered (thus polarized) by electrons and dust into our line of sight. 

The nucleus of the galaxy is not well defined at UV and visible wavelengths; instead, high resolution 
UV and visible imaging with the Hubble Space Telescope (HST) only shows a number of bright clouds and 
filaments (Evans et al. 1991; Macchetto et al. 1994; Bruhweiler et al. 2001). There have been a number 
of attempts to estimate its location from the direction of the UV polarization vectors measured with HST 
(e.g., Capetti et al. 1995a, 1995b; Kishimoto 1999). On the other hand, near-infrared (NIR) and mid- 
infrared images (Braatz et al. 1993; Marco, AUoin, & Beuzit 1997; Thatte et al. 1997; Rouan et al. 1998; 
Weinberger, Neugebauer, & Matthews 1999; Bock et al. 2000) show an intense unresolved source. Thompson 
et al. (2001), by correlating the optical Ha and NIR Pa images from HST, determined that the NIR point 
source lies at the location of the nucleus (to within the errors of measurement) as determined by Kishimoto 
(1999). These locations can be seen in Figure 7 and Table 2 of Kishimoto (1999), where he plots the nucleus 
position determined from the UV polarization measurements on the HST FbOlN [O III] image. The position 
(with imcertainty of ^ 0.04") of Thompson et al. is consistent with both Kishimoto's nucleus position and 
with Cloud B of Evans et al. (1991), also plotted on the figure. The radio source SI is thought to be the 
location of the AGN nucleus and accretion disk because of the high velocities of its H2O masers emanating 
from an apparently rotating, 4 pc diameter disk (Gallimore et al. 1996b, 2001; Greenhill & Gwinn 1997) 



Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with GTO proposal #7212, PI: E. F. Erickson. 
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with interior hot thermal emission (Muxlow ct al. 1996; Roy et al. 1998). The coordinates of source SI 
agree with the positions of Cloud B and the nuclear position plotted by Kishimoto (1999) within the absolute 
uncertainty of 0.08" of the optical position of the galaxy (Capetti, Macchetto, & Lattanzi 1997b). 

If the intense NIR unresolved source is really the galaxy and AGN nucleus, its light should be scat- 
tered and polarized with a centrosymmetric polarization pattern. Moreover, since NIR wavelengths are less 
affected by extinction than visible or UV wavelengths, the morphology of the scattering region can be bet- 
ter determined through measurement of the NIR polarized light. Ground-based measurements of the NIR 
polarization in > 1" to 2" resolution have been made by Young et al. (1996), Packham et al. (1997), and 
Lumsden et al. (1999). They show that the polarization is not consistent with scattering in the immediate 
vicinity of the nucleus; Young ct al., Packham ct al., and Lumsden et al. all suggest that there is additional 
polarization caused by the light from the nucleus being absorbed by aligned dust grains in the AGN torus. 
Young et al. and Lumsden et al. also detect a dip in the polarized emission ~ 1" south of the nucleus and 
all three groups see a southwest lobe of polarized emission that is not observed at visible wavelengths. 

In this paper we report on 2 /xm polarization observations taken by the Near Infrared Camera and 
Multi-Object Spectrometer (NICMOS) on HST. With the 0.2" spatial resolution of NICMOS, we sec the 
unresolved polarized source that must be the nucleus and numerous scattering clouds in the NE and SW 
lobes. Section 2 describes the observations, section 3 the results and discussion, and section 4 summarizes 
the results. 



2. Observations and Data Reduction 

NGC 1068 was observed in two visits on 11 October and 15 November, 1997, with the telescope rolled 
such that NICMOS Camera 2 was oriented at 63.519° and 154.519° East of North, respectively. Integrations 
were taken at 4 dither positions separated by 10.5 pixels (0.78") in each of the three Camera 2 polarizing 
filters, POLOL, POL120L, and POL240L, which cover a 1.9 to 2.1 /im bandpass. The detector array was 
read out in MULTIACCUM mode with sample sequence STEP16 to accumulate total times of 64 or 80 
seconds per integration for visit 1 and 128 seconds for visit 2. We were able to recover the full dynamic 
range in this observation even though the nucleus saturates in 8 or 16 s because the design of the STEP16 
MULTIACCUM sequence incorporates very short dwell times for the first few read-outs. 

The raw data were calibrated using CALNICA (v. 3.3) as implemented by STScI in STSDAS v2.2. We 
included the CALNICA task BIASEQ but not the tasks PEDSKY or PEDSUB because there is no blank sky 
in the 19.4" images. The reference files used were the most recent Institute flats and non-linearity files (of 
January, 2001), but the dark files were provided by M. Rieke (2000, private communication). Each reduced 
image required additional corrections for quadrant bias and half the images needed slight adjustments for 
residual shading. Since images requiring no shading adjustment occurred during both visits, using the STScI 
temperature sensitive darks would not have negated the need for adjustments. These corrections were very 
small, < 0.05 counts whereas the nucleus center is ~ 3600 counts s~-^ and the image edges ~ 0.5 counts 
s^^. Although the need for such corrections is not noticeable in the total intensity images, it is obvious in 
quantities that incorporate a difference image, such as the polarized flux. 

For each polarizing filter, we aligned and shifted the four dither positions by centroiding the unresolved 
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nucleus and then median combined them to remove bad pixels using the IDL program, IDP3^. We used the 
equations of Hines, Schmidt, & Schneider (2000) to calculate the total flux (F), percentage polarization {p), 
polarized flux {pxF), and position angle for each visit; we calculated statistical errors using the equations of 
Sparks & Axon (1999) and the estimated errors on the flux computed by CALNICA. The original pixel scale 
was 0.0760" by 0.0753" per pixel; the pixels were rectified and smoothed with a 3 x 3 smoothing function 
(i.e., to approximately the ffST spatial resolution of 0.2" at 2.0 /um) in the data reduction process to achieve 
higher signal/noise. 

The total intensity, fractional polarization, and polarized flux images from the two visits are plotted 
in Figure 1. Two ghosts at pixel positions (23,35) and (8,20) were discovered in the POLOL and POL120L 
filters, respectively, and are marked in the images. Although the total fluxes for the two visits arc very 
similar, the percentage polarization and polarized flux differ by 2-3 percent with visit 2 generally showing 
less polarization. The cause of the difference in polarization is not known. Because the differences between 
the visits in the percentage polarization and polarized flux are generally larger than the calculated statistical 
errors, we estimated the actual uncertainties (which we quote in later sections) as equal to 0.5 times the 
differences between the visits. 



3. Results and Discussion 
3.1. Nucleus 

In Figure 1 the nucleus of NGC 1068 appears as an intense, unresolved source convolved with the 
HST point spread function (PSF) in both the total flux and the polarized flux images. Since we do not 
have a real-time example of the PSF from a star, we calculated the PSF for each of the visits using the 
program TinyTIM^ (Krist & Hook 1997, 1999). The nuclear region with contours of the Tiny TIM PSF is 
plotted in Figure 2. It is not possible to remove the effects of the PSF by subtracting the calculated PSF 
from the observed nucleus — numerous blobs at the locations of PSF diffraction rings and spikes remain 
no matter what PSF scaling was used. In fact, one can sec in Figure 2 that although there is general 
agreement between the calculated PSF contours and the polarized flux, there are also many regions that do 
not have good enough correspondence between the calculated PSF and the observed blobs for complete PSF 
subtraction. As a result, we can only say that there is a hint of additional structure within the first Airy ring 
surrounding the nucleus. It is possible that the nucleus is slightly extended in the north-south direction as 
was suggested by Weinberger et al. (1999), Rouan et al. (1998), and Thatte et al. (1997). It is interesting 
that apparently the radio jet is initially emitted in approximately a north-south direction and is diverted to 
the ~ 30° position angle 0.3" north of the nuclear source SI (Gallimore, Baum, & O'Dea 1996a; Muxlow et 
al. 1996). 

The nuclear polarization, averaged over the central 3x3 pixels, is 6.0 ± 1.2% and the position angle 
is 122° ± 1.5°, perpendicular to the eventual direction of the radio jet but not its initial direction. This 
percentage polarization is larger than that measured by Packham et al. (1997) or Lumsden et al (1999), 
probably because their measurements were diluted by the non-polarized flux from the galaxy core seen in 
their larger beams. The measured position angles are very similar. 



^http: / /nicmos. as.arizona.edu 

*http: / /www.stsci.edu/softwcire/tinytim 
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Lumsden et al. (1997) studied their polarization position angles to check for the centrosymmetric pattern 
that would indicate scattering of light originating at a well defined location (e.g., the nucleus). They found 
that the NIR polarization vectors in the vicinity of the nucleus are not consistent with scattering of nuclear 
light only; instead they suggested that the nuclear emission is polarized by dichroic absorption as a result 
of passing through aligned dust grains. However, their beams were typically 2" or larger and thus could 
include non-nuclear light as well as nuclear. Alternatively, the nuclear polarized light could include light 
backscattered from the cones above the nucleus and scattered a second time on the accretion disk (Fischer, 
Henning, & Yorke 1996; Schinnerer et al. 2000). 

In order to estimate what fraction of the nuclear polarized light is coming from the point source, we 
measured the nuclear fluxes and the calculated PSF in various apertures and corrected the fluxes by dividing 
by the PSF fraction in the same apertures. These results are given in Table 1. We see that the polarized flux 
from the nucleus of NGC 1068 completely dominates the polarized flux that would be seen in any 1" — 2" 
aperture such as were used in the ground-based polarization measurements; contributions from off-nuclear 
scattered light would be barely detectable. Only a small fraction of the polarized flux measured within a 
2" aperture is not flux from the essentially unresolved nucleus (subtracting the PSF-corrected polarized flux 
measured in the 0.209" radius aperture from that measured in the 1.00" radius aperture gives 3.4±0.2 mJy). 

The polarization vectors are plotted on a contour map of the polarized flux in Figure 3. Since the vectors 
describing the vibration of the electric field in the sky plane are perpendicular to the projected fine between 

the illuminating source and the (optically thin) scattercrs, one can locate the illuminator by computing 
the intersection of the lines perpendicular to the polarization position angle vectors (e.g., Lumsden et al. 
1999; Kishimoto 1999; Capetti et al. 1995a). The procedure is to pick a pixel, assume it is the center, 
calculate the vector angle to every other pixel, and subtract the perpendicular polarization angle at those 
pixels. The means of the differences as a function of position angle from the assumed center are a minimum 
when the assumed center pixel is located at the illuminating source of the scattered light. For the NICMOS 
polarization vectors for regions greater than 1" from the peak, the polarization vectors point to the intense 
peak within < 1 pixel. This confirms that the illuminating source of the scattered light is the same at 2 fim 
as the illuminating source estimated from optical polarization measurements (e.g., Kishimoto 1999), given 
the alignment of the optical and Paschen-a images of Thompson et al. (2001). 

3.2. Inner Two Arcseconds 

Even though most of the brighter emission blobs in the central 2" arc due to the complex diffraction 
pattern from the bright central point source (i.e., the PSF), it is clear that there is some polarized emission 
extending along a ridge < 1" north and south of the nucleus (Figs. 2 and 4). This appears to be the structure 
imaged by Thompson et al. (2001) in Pa (their Fig. Ic) and may include some of cloud C of Evans et al. 
(1991). The other clouds seen in the inner 1" in UV polarization (e.g., Capetti et al. 1995a, 1995b) are 
obscured by the PSF and there is no good correspondence with the other UV polarization structures in the 
central 2". The emission spot described by Thompson et al. (2001) is not bright in polarized light. 

Young et al. (1996) and Lumsden et al. (1999) noticed that there is a dip in the polarized flux about 
1" south of the nucleus in their J and H images but not in their K band image. This dip is clearly present 
in our 2 /xm NICMOS image (see especially the fractional polarization images in Fig. 1), although it is not 
clear whether it is a completely empty gap or whether there are some clouds of scatterers present (the fact 
that the observed polarized flux blobs do not completely correspond to the Tiny TIM- calculated PSF is to be 
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expected and does not necessarily indicate additional scattering sources in the region). Young et al. (1996) 
suggested that the dip is caused by obscuration from the dusty torus surrounding the AGN; this torus would 
have to have a size ^ 200 pc, which is substantially larger than the few-parsec-sized canonical AGN torus 
(see Antonucci 1993 and references therein) that directs the radiation into the observed scattering cones. 
Another possible source of the obscuring dust is the dense cloud or nuclear ring (diameter ~ 5") seen in CO 
by Schinnerer et al. (2000). Schinnerer et al. estimate this ring to have a mass of 5 x 10^ M© centered on 
the nucleus. Using the interstellar ratio of Ay/hydrogen column density of Bohlin, Savage, & Drake (1978), 
we estimate values of Ay for this cloud ranging from 6 to 15. Although this cloud could provide the small 
dip seen by Young et al. and Lumsden et al., it would not provide enough extinction to produce a gap as 
empty as the polarization gap in Fig. 1 appears to be. We suggest that there is a real dearth of scatterers 
in the band 1" south of the nucleus. 



3.3. Southwest and Northeast Lobes 

The southwest and northeast lobes (sec Fig. 4 for nomenclature) are seen to be highly clumped, with 
locally very strong polarization of 17% when smoothed with a 3 x 3 smoothing to 0.23" per smoothed pixel 
(the NICMOS resolution is 0.2") and 20% in a single pixel. The southwest lobe is quite faint at UV and 
visible wavelengths; the reason probably is that it lies behind the plane of the host galaxy and is extinguished 
by interstellar dust. Observations of the Ha/H/3 ratio by Bruhweiler et al. (2001) show that the northeast 
lobe seen in the ionized lines suffers very little extinction but the southwest lobe has Ha/H/? > 6 (the Case 
B recombination line ratio is 2.87). This implies Ay > 2 mag for ionized hydrogen clouds to the southwest 
of the nucleus. 

In addition to inferring the location of the main illuminating source, one can also use the polarization 
vectors to infer whether there is any additional illuminating source besides the strong intensity peak (which 
is very close to the illuminating source, if not the actual source), or whether there is additional polarization 
from dichroic absorption by aligned interstellar dust grains in the galaxy. To do this, we measured the 
deviations of the perpendiculars to the polarization vectors from the vectors to the bright nucleus for regions 
not containing PSF contributions from the nucleus. Most of the largest deviations from zero occur where 
the intensity is less than the local maximum, in other words, in the regions of larger uncertainty. The only 
locations with large deviations that are consistent between the visits are at the extreme west edges of the 
northern and southern lobes. The uncertainty, though, is of similar size to the deviation, about 2 degrees 
whereas the deviation is about 4 degrees. We conclude that there is no strong evidence of any additional 
polarization component in the regions 1" — 4" from the nucleus from dichroic absorption of grains aligned 
by a magnetic field oriented other than perpendicular to the line joining the nucleus and the grain location 
in the galaxy. However, since magnetic fields in galaxies, including NGC 1068 (Scarrott et al. 1991), tend 
to follow the spiral arms, it is likely that any dichroic absorption would produce polarization with the same 
orientation as that produced by scattering of light emitted by the nucleus. On the other hand, we find 
that there is no correlation of the bands of polarized flux in the southwest lobe and regions of high Ha/H/3 
(large extinction) in Figure 3 of Bruhweiler et al. (2001). Thus, dichroic absorption is probably not the 
main contributor to the high polarization that we see in the southwest lobe in spite of the high optical 
extinction with the possible exception of the extreme southwest clump, where Bruhweiler et al. (2001) saw 
no appreciable Ha emission. Most likely the higher polarized intensity in the southwest cone compared to 
the northeast cone is that it contains a higher density of dust (larger filling factor). 

Jets emitting synchrotron radiation have been observed emanating from the nucleus by Wilson & Ulves- 
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tad (1983, 1987). In Fig. 4, wc plot the 4.9 GHz radio contours on top of our combined polarized flux image. 
The jets appear to pass through the regions containing the scattering material without spreading out; only 
after they are beyond the regions exhibiting polarized emission do they spread out into extended lobes. This 
is the first observed correlation of some aspect of the south radio lobe with any optical, IR, or X-ray feature. 
One possibility is that the scattering material in the NE and SW lobes seen in NIR polarization is affecting 
the radio jets and the formation of the radio lobes; another is that the jets have pushed material into the 
dense dust clouds that are visible in scattered light. 

3.4. Northecist Knot 

The northeast "knot" was first identified by Miller et al. (1991) as a region or knot of high visible 
polarization about 4.7" northeast of the nucleus. An enhancement in the polarized emission is also apparent 
in the NIR images of Young et al. (1996), Packham et al. (1997), and Lumsden et al. (1999). However, 
this region does not stand out in the optical ionized line images of Bruhweiler et al. (2001): ionized gas is 
present but extended over a much larger region than the northeast knot, the excitation from the [O III]/II/3 
ratio is much lower right at the position of the knot than the surrounding gas, and the extinction inferred 
from the Ha/H/? ratio is only moderate. 

The northeast knot is intriguing because it is clearly on the outside edge of the NE radio lobe observed 
by Wilson & Ulvestad (1983, 1987). In fact, it looks like the jet is forced to go past the cloud/knot and 
cannot penetrate it, as though the cloud is a massive obstruction in the path of the radio lobe. There is a 
fainter cloud of polarized emission immediately north of the radio lobe as well. Wilson & Ulvestad (1987) 
suggested that the radio polarization morphology and limb brightening of the NE radio lobe are caused by 
the high energy particles of the lobe encountering the interstellar medium of the galaxy in a radiative bow 
shock. Schinnerer et al. (2000) measured ^^CO (1 0) emission about 5" north of the nucleus. Their largest 
cloud is just east of the location of the northeast knot — we suggest that the northeast knot and the north 
cloud are the edges of this cloud that is facing the AGN nucleus. This, then, would be the cloud that is 
confining the NE radio lobe. Miller et al. (1991) also suggested that the northeast knot has a direct view of 
the nucleus and its ionizing radiation. 

3.5. Electron Scattering vs. Dust Scattering 

It has been thought for some time that the ultraviolet and visible polarized flux in the nucleus of 
NGC 1068 is due to scattering of light by electrons in the vicinity of the AGN central engine (Miller et al. 

1991; Antonucci 1993; Yoimg et al. 1995). Even though scattering by dust is more efficient than electron 
scattering, electron scattering is suggested for the nuclear region because of the constancy of percentage 
polarization between UV and visible wavelengths, because the X-rays appear to be scattered in the same 
way as the visible polarized flux (Antonucci 1993), and because the H/3 line shows evidence of broadening 
from scattering by high temperature electrons (Miller et al. 1991). The Chandra X-ray images (Young, 
Wilson, & Shopbell 2001) show good correlation with [O III] and Ha images (Capetti, Axon, & Macchetto 
1997a; Macchetto et al. 1994; Bruhweiler et al. 2001), which are similar to the HSTUV polarized intensity 
images (Capetti et al. 1995a, 1995b). 

If the polarized flux from the unresolved nucleus is due to dichroic absorption, as seems likely, then 
we should compare the non-nuclear NIR polarized flux to the visible polarized flux to investigate whether 
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both are due to electron scattering. Inglis et al. (1995) measured the polarized spectrum from 0.43 to 0.69 
/xm in a 1.3" x 2.7" slit at a number of positions, including positions 2.5" and 5.0" NE of the nucleus. In 
the position 2.5" NE of the nucleus (approximately the NE lobe), their polarized flux decreased from about 
0.4 mJy to about 0.2 mJy over this spectral range. Our measured polarized flux in the same aperture is 
> 0.54 mJy after wc subtract 0.15 mJy for the contribution from the diffraction spikes from the nucleus, 
as computed using Tiny TIM (the actual peak of the NE lobe is a little further to the West and has larger 
flux). The position 5.0" NE of the nucleus is that of the northeast knot. Here Inglis ct al. (1995) found that 
the polarized flux decreased from about 0.35 mJy at 0.45 /xm to about 0.25 mJy at 0.69 fim. In a similar 
aperture we measured the polarized flux to be 0.65 mJy in the northeast knot at 2.0 /im. For both positions 
the ratio of the polarized flux at 2.0 /xm divided by the polarized flux at 0.6 fim is the same, ~ 2.2. We 
conclude that the scattering is produced by the same process at both positions. 

The possibilities for the large increase in scattered flux from visible to NIR wavelengths at these two 
non-nuclear positions are: 1) The intrinsic NGC 1068 Seyfert 1 AGN spectrum increases rapidly between 
visible and NIR wavelengths, even more so than is typical (Ward et al. 1987 observed that the fluxes, 
of bare, unreddened Seyfert 1 galaxy nuclei increase from visible to NIR wavelengths even though they 
are flat or increase from visible to UV wavelengths). Such an increase in the nuclear flux could be due 
to an increasing contribution from the very hot dust immediately surrounding the central engine. 2) The 
polarized flux seen at visible wavelengths suffers substantial extinction such that the apparent flux ratios of 
NIR/visible are much larger than the intrinsic ratios. If this were the case, though, there should be almost 
no UV polarized flux and the slowly decreasing polarized flux as a function of visible wavelength (Miller et 
al. 1991; Inglis et al. 1995) would not be observed. 3) There are systematic effects such that NICMOS 
measures too much polarized flux. Overall, however, our nuclear photometry at 2.0 fim is consistent with 
other NIR measurements (Weinberger et al. 1999; Lumsden et al. 1999; Packham et al. 1997), given the 
steep rise in nuclear flux between the H and K bands, and thus we do not consider this possibility to be 
likely. 

For the northeast knot. Miller et al. (1991) suggested that the polarization was caused by dust scattering 
because the percentage polarization decreases rapidly to long wavelengths over the observed wavelength range 
of 0.33 to 0.65 /im and because the polarized H/9 line is narrower than the nuclear polarized H/3 line. An 
additional reason for rejecting electron scattering as the cause of the polarization in the northeast knot is that 

the northeast knot is not particularly prominent in the Chandra image (Young ct al. 2001), whereas electron 
scattering is suggested as dominating the nuclear polarized flux because of the similarity of appearance of 
the nucleus in UV and X-ray images. If the scattering is also caused by dust at 2 /xm, the nuclear source flux 
would have to increase by much more than a factor of 2.2 since the scattering cross section for interstellar 
dust decreases by a factor of order 20 from 0.6 /im to 2.0 /im. Inglis ct al. (1995) modeled the scattering and 
polarization in the central regions of NGC 1068 with a model that includes both electron and dust scattering 
(Young et al. 1995). They could not get a fit to their NE positions unless they included electron scatterers 
as well as dust in the dust cloud at the location of the northeast knot. Packham et al. (1997) also modeled 
the scattered emission from the northeast knot including both dust and electron scattering and found an 
acceptable fit. Electrons are certainly present, given the ionization implied by the HST ionized line images 
(e.g., Capetti et al. 1997a; Bruhweiler et al. 2001; Thompson et al. 2001) and optical spectropolarimetry 
(Miller et al. 1991; Inglis et al. 1995). This ionization is from photons from the AGN and not from shocks 
from the radio lobes (Bruhweiler et al. 2001). It is probable that both electrons and dust contribute to the 
NIR scattered light from these non-nuclear regions. 
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4. Summary 

We have used HST NICMOS Camera 2 to observe the polarized Hght in the center of NGC 1068 at 
2 iim with a difFraction-hmited PSF of 0.2" FWHM. These results are one more example of the value of 
high-spatial-rcsolution observations in disentangling complex sources, here the nucleus of one of the most 
important AGNs. Since the clouds emitting the polarized light measured in this study see the AGN directly, 
future work with a similar high-spatial-resolution instrument might be to observe any nuclear activity through 
its echos in the polarized flux as a function of nuclear distance (light travel time) (e.g., Gallimore et al. 2001). 

The nucleus is an intense unresolved source, polarized at a level of 6% with a position angle of 122°. It 
is likely that the polarization is due to dichroic absorption in the overlying molecular cloud. 

There are two polarized lobes extending up to 8" northeast and southwest of the nucleus. The polarized 
flux in both lobes is very clumpy, with the maximum polarization occurring in the southwest lobe at a level 
of 17% when smoothed to 0.23" resolution. The southwest lobe is probably located behind the plane of the 
galaxy because it cannot be detected at optical wavelengths in spite of the fact that it is much brighter than 
the northeast lobe at 2 fiia in polarized flux. The perpendiculars to the polarization vectors in these two 
lobes all point back to the intense unresolved source within one 0.076" Camera 2 pixel, thereby confirming 
that this is the illuminating source of the scattered light and therefore the probable AGN central engine. 

Features in the polarized lobes include a gap at about 1" between the nucleus and the southwest lobe 
and a "knot" of emission about 4.5" northeast of the nucleus. Both features had been discussed by previous 
authors, but they are much better defined with the high spatial resolution of NICMOS. The gap is sufficiently 
deep that we suggest that it is due to a lack of scatterers rather than absorption by the overlying molecular 
cloud seen in CO by Schinnerer et al. (2000). The polarized flux images were aligned with the radio map 
of Wilson & Ulvestad (1983) by assuming that the intense point source is located at the same coordinates 
as the radio source SI (it is within 0.1" of SI according to the measurements of Capetti et al. 1997b and 
Thompson et al. 2001). The northeast knot and the north cloud may be the side of a molecular cloud that 
is facing the nucleus. This cloud may be constraining the northeast radio lobe at the head of the radio 
synchrotron-emitting jet. 

Polarized fluxes were measured for positions 2.5" and 5.0" NE of the nucleus in a 1.3 x 2.7" aperture 
to correspond to the positions measured by Inglis et al. (1995) from 0.43 to 0.69 /im. (The nucleus position 
was not considered because of the severe effects of extinction at the shorter wavelengths and excess polarized 
flux from dichroic absorption at 2 /um.) The ratio of 2 /xm flux to 0.6 /xm flux is ~ 2.2 at both positions. 
This could indicate that the same scattering processes obtain at both 2.5" and 5.0" from the nucleus. It 
also shows that the intrinsic spectrum of the NGC 1068 nucleus increases strongly from optical wavelengths 
to NIR wavelengths, making it difficult to determine the wavelength dependence of the cross sections of the 
materials doing the scattering, since the intrinsic spectrum of the source of the scattered light — the AGN 
central engine — is not known. 
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Fig. 1 . — The logarithm of the total intensity, the fractional polarization, and the logarithm of the polarized 
intensity at 2 fim are plotted for both visits. The intensities are scaled so that the nucleus is saturated. 
The large diamond is the POLOL ghost and the large triangle is the POL120L ghost. Figure la-c: Visit 1. 
Figure Id-f: Visit 2. The Visit 2 images have been rotated by 90° so that they have approximately the same 
orientation as the Visit 1 images. 
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Relative Right Ascension (arcsec) Relative Right Ascension (arcsec) 

Fig. 2. — The unsmoothed log polarized flux images (same flux scale as Fig. 1) from the two visits are 
overlaid with contours from the TinyTIM PSF calculation. The contour levels are at 0.00032 (dot-dashed), 
0.001, 0.0032, 0.01, 0.032, 0.1, and 0.32 of the peak. The ghost in the POL120L filter is marked by a triangle. 
Both images have been rotated so that North is Up. Features of interest that are not due to the PSF include 
the bright unresolved nucleus (magenta), the ridge of polarized flux (red) extending about 1" due north of 
the nucleus, the gap of low polarized flux (green and yellow) about 1" south of the nucleus, and the other 
regions of polarized flux (yellow and orange) that lie outside of any PSF contour, especially the bands of 
polarized flux about 2" south of the nucleus in the SW Lobe. The areas that appear to be polarized flux 
WNW and ESE of the nucleus are the i75T diffraction spikes from the nucleus, a: Visit 1. b: Visit 2. 
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6 4 2 -2 -4 -6 
Relative Right Ascension (arcsec) 

Fig. 3. — Contours of the polarized intensity, smoothed by 5 pixels (0.38"), are plotted, along with the 
polarization vectors. North is up. The contour levels are 0.08, 0.16, 0.32, 0.64, 1.28, 2.56, 5.12, 10.24, 20.48, 
and 40.96 mJy/square arcsec. The unmarked maximum intensity is 76.3 mJy/square arcsec. The large 
diamonds mark the POLOL ghost and the large triangles mark the POL120L ghost. 
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Fig. 4. — The combined polarized flux images from the two visits are overlaid with contours of the 4.9 GHz 
radio measurements of Wilson & Ulvestad (1983). North is up. The radio contour levels are 0.0005, 0.001, 
0.002, 0.004, 0.008, 0.032, 0.128, and 0.256 Jy beam"^ The images are registered by assuming that the 
source SI (Muxlow et al. 1996; Gallimore et al. 1996b) is the AGN nucleus and is at the location of the 
intense NIR peak. At the resolution of the radio image (0.4" x 0.4"), the radio peak is dominated by radio 
source C, 0.29" north of source SI (Roy et al. 1998), and source SI is not apparent. 
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Table 1. Photometry of NGC 1068. 



Location 


Aperture 


PSF* 


Measured 


PSF-Corrcctcd 


Polarized 


PSF-Corrected 




arcsec 




Flux (Jy) 


Flux (Jy) 


Flux(mJy) 


Polarized Flux (mJy) 


Nucleus 


0.209^ 


0.531 


0.188 ibO.OOr 


0.353 ±0.002 


11±2 


20 ±4 


Nucleus 




0.831 


0.315 ±0.001 


0.379 ±0.001 


18 ±4 


21 ±4 


Nucleus 


1.000<i 


0.908 


0.388 ± 0.001 


0.427 ±0.002 


22 ±4 


24 ±4 


NE Knot 


1.9 X 1.9** 








0.74 ±0.13 





"^Fraction of PSF within the aperture calculated by TinyTIM 
''Radius of circular aperture corresponding to the first AIRY dark ring 
'^Errors calculated by dividing the difference between the visits by 2 
''Radius of circular aperture 
®Size of square aperture 



